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ABSTRACT 



PUMPING-EROSION SUBSIDENCE 
OF A SEAFLOOR PLATE-FOOTING 
(November 1987) 

ALVIN EUGENE GRIMMIG JR., B.S. University of Washington 
Chairman of Advisory Committee: Dr. John B. Herbich 

Subsidence of objects that rest on the seafloor is an 
expected but not well understood phenomenon. Subsidence is 
the settlement of fixed structures and mobile structures 
(such as oceanographic packages) that are placed on the 
seafloor. A major concern involving bottom resting systems, 
in relatively shallow water (less than 200 meters) is wave 
induced motions. The system response of the structure in 
this dynamic environment can create substantial loadings at 
the structure-seafloor interface. Past works have addressed 
wave-induced and current scour as the major factors causing 
structural subsidence, neglecting the effects of pumping - 
erosion. 

The research described is directed towards indentifying 
and understanding the phenomenological aspects of pumping- 
erosion. The work describes a new hypothesis concerning 
liquefaction as being the precursor to subsidence in a two 
cycle pumping-erosion process. Due to the complexity of 



structure-soil interaction a multivariable linear regression 
was performed on the data to substantiate experimental 
observations. 



ACKNOWLEDGEMENTS 



The author is grateful to Dr. John B. Herbich, for 
serving as his Committee Chairman. Appreciation is also 
extended to Dr. Wayne Dunlap and Dr. Wilfred Gardner for 
their role as Committee Members and for their advice, 
interest and encouragement. A special thanks to Dr. Robert 
Randall, for his perspective and encouragement. 

I am especially grateful to wife, Debra and sons Matt, 
Andy and Aaron, for their continuous love, faith and support 
during this course of study. In appreciation, I dedicate 
this work to them. 

Thanks is also extended to technicians Randy Bush, 
Charles Carnes, and Carl Fredrickson for their assistance in 
the development and operation of equipment during research. 



TABLE OF CONTENTS 



ABSTRACT iii 

ACKNOWLEDGEMENTS iv 

TABLE OF CONTENTS v 

LIST OF TABLES vi 

LIST OF FIGURES vi i 

CHAPTER I INTRODUCTION 1 

CHAPTER II LITERATURE REVIEW 7 

CHAPTER III PROCEDURES AND INSTRUMENTATION 11 

CHAPTER IV RESULTS 19 

CHAPTER V DISCUSSION OF RESULTS 34 

CHAPTER VI SUMMARY and CONCLUSIONS 65 

REFERENCES 69 

APPENDIX A - SUBSIDENCE CURVES 71 

APPENDIX B - EXPERIMENTAL RECORD (SAMPLE) 104 

APPENDIX C - SCOUR CONTOURS 110 

VITA 154 



LIST OF TABLES 



Table Page 

1. Experimental Parameters 20 

2. Submerged Angle of Repose 43 

3. Variable Range of Values 48 

4. Variable Correlations 49 

5. Order of Variable Significance 52 

6. PROC REG - 8 variable 55 

7. PROC REG - 6 variable 56 

8. PROC REG - combined factors 57 

9. Multivariable Model Summary 63 



LIST OF FIGURES 



Figure Page 

1. Structure-Soil Response of a Seafloor Plate-Footing.. 3 

2. Sediment Test Tank 12 

3. Grain Size Analysis 13 

4. Experimental Set-Up 15 

5 ( A-L ) . Typical Strip Chart Record of Pore Water Pressures 
and Exit Velocities 22 

6. Upward Plate Cycle 36 

7. Distorted Grain Skeleton 37 

8. Downward Plate Cycle 39 

9. Bed Profile Details 40 

10. Post Negative Pore Water Spikes 42 

11. Ranking of Experiments by Liquefacton Index 45 

12. Ranking of Experiments by ( A2Ap/WT ) 46 

13. Predicted (8-variable) Model vs Actual 60 

14. Predicted (6-variable) Model vs Actual 61 

15. Predicted (2-variable) Model vs Actual 



62 



CHAPTER I 



INTRODUCTION 

Structural subsidence can range from a few inches to 
several feet causing results varying from slight tilting and 
displacements of a structure to total and catastrophic 
structural failure. As a result of the petroleum and gas 
industries use of seafloor supported structures (i.e. jack- 
up rigs) much investigation has been directed at the 
structure-soil system and interaction. Although this area 
has been intensely investigated, the focus of attention has 
been on large oil rig structures. 

The use of relatively small mobile submerged structures 
for oceanographic and industrial purposes has increased in 
the past several years. The ability to successfully deploy, 
operate and retrieve these instruments and equipment 
packages is largely dependent upon the structure-seafloor 
interaction. Bottom resting systems are subject to the 
dynamic response of the structure in the ocean regime which 
can produce significant dynamic loadings at the seafloor. 

The stability of a submerged seafloor structure is dependent 
upon the subsidence and orientation of the structure while 
subject to wave and current forces. These forces produce a 
footing response in a variety of bearing, tension, sliding. 



rocking and torsional strains in the marine soil supporting 
the structure (Figure 1). 

The central issue of this report involves the 
investigation of the phenomenon and mechanisms of pumping- 
erosion and its contribution to subsidence. Pumping-erosion 
is inherently non-linear as the structure is intermittently 
coupled and decoupled from the soil interface. In addition, 
the effects of the period of oscillation, amplitude of 
displacement, footing characteristics, hydrodynamic effects 
and soil parameters are simultaneously involved resulting in 
a very complicated and complex problem. 

The basic pumping-erosion cycle involves an uplifting 
force vertically displacing the footing off the seafloor. 
Subsequently, a downward force will dynamically load the 
structure onto the soil again. This over simplified cycle 
is a Trojan horse, as it involves several discrete and 
intertwined forces which are not readily apparent. 

Consideration must be given to the geometry and 
dynamics of the footing, The size, shape, and weight of the 
structure will determine such characteristics as center of 
mass, center of buoyancy, drag and lift forces and moments 
of inertia which will determine the structures response to 
the exciting force. The wave and current environment will 
produce the exciting force characteri sties of period. 
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amplitude, and duration of excitation. Hydrodynamical ly , 
effects of exit and entrance velocities, water density and 
viscosity as well as the added mass of the interacting 
structure will influence the structure response and the 
sediment transport and settlement. 

Soil parameters such as soil density, permeability, 
void ratio, grain size and shape, and pore water pressures 
will determine seafloor response to the structure-soil 
interaction. 

Traditionally subsidence was attributed to the effects 
of scour resulting from the driving force of wave and 
currents with little regard to pumping-erosion. Any 
inference to pumping-erosion usually surrounds the latter 
half cycle where the water is displaced or "pumped" from 
under the footing as it makes contact with the soil surface. 
The soil transport has been attributed to the "jetting" 
action of the water escaping from under the footing. 

But what of the first half of the cycle of the footing 
when it is uplifted from the seafloor ? What forces are 
involved ? And what contributions do these forces make to 
the process of subsidence ? As the footing is uplifted, 
negative pore water pressures can be developed in the soil 
under the footing, possibly to the point of liquefaction. 



Liquefaction, or fluidization, occurs when the pore- 
water pressure in a soil is equivalent to the total stress 
of the soil; in other words when the effective strength of 
the soil equals zero. Liquefaction of a soil requires 
considerable pore water pressures as compared to the 
relative ease of fluidization, which reduces the effective 
strength of the soil by the increase in the pore water 
pressure combined with the uplifting effect of the upward 
seepage flow. If the bed is considered deformable, the 
fluidization would then weaken the grain skeleton and 
distort the grain bed surface causing the grains to have a 
greater hydraulic profile making them readily available for 
transport during the downward half cycle of the footing. 

Thus, it is hypothesized that subsidence is a two phase 
cycle. The first phase occurs on the uplifting of the 
footing causing sufficient negative pore water pressures to 
fluidize the bed which weakens and distorts the bed making 
the upper grain layers more susceptible to sediment 
transport. The second phase, occurs on the downward stroke 
of the footing where the velocities developed under the 
footing transport the sediment. The transport to the 
sediment will be accentuated at the edges where the water is 
jetted out from under the footing. The following research 
is designed to investigate the hypothesis on a 
phenomenological basis. 



OBJECTIVES 



The objectives of this research are as follows: 

1. To determine which geometric, 
hydrodynamic and kinematic parameters are involved 
in pumping-erosion subsidence. 

2. To evaluate the relative importance of 
the geometric, fluid and kinematic parameters 
involved in pumping-erosion subsidence. 

3. To determine if sufficient negative pore 
water pressures can be developed under the footing 
to fluidize the foundation bed to the point of 

1 iquefaction. 



CHAPTER II 



LITERATURE REVIEW 



There are a limited number of references dealing with 
soil/water/structure interaction. Even fewer address the 
soil /water/structure interactions and forces involved in 
structural subsidence. Furthermore, there are only two 
sited works that discuss the topic of pumping-erosion. Most 
structural subsidence is attributed to he effects of scour. 
Herbich et a 1 . ( 3 ) has compiled an extensive biography on the 
topic of conventional scour. 

Reimnitz and Kempema (10) in 1981, describes the 
effects of dynamic ice wallowing in the Alaskan nearshore 
areas. The report hypothesizes the occurence oftwo types of 
hydraulic processes. In the first process the ice plays a 
passive role, acting as a flow obstacle in the current and 
wave regime. This particular aspect would account for the 
more conventional scour around structures. The second 
hydraulic process involves ice playing an active role, 
either by simple vertical oscillations or by wallowing 
(rocking). Dahlberg (1) discusses pumping erosion as it 
pertains to structures. He states that pumping erosion is 
associated with excess pore-water pressures set-up in the 
foundation soil during storm periods. Dalhberg suggests 



that excess pore-water pressures develop to balance the 
overturning moments on the base of the structure. It is 
further stated that if there is free communication between 
the foundation soil and the seabed, the pore-water gradient 
may be high enough to liquefy the soil locally. 

It was verified that a North Sea platform, Frigg CDPI, 
which had no protective scour skirts, was subject to 
pumping-erosion. Divers observed periodic puffs of sediment 
around the periphery of the structure base when the 
sediments were carried in suspension as the structure rocked 
in the storm environment. It is noted that the pumping- 
erosion described was solely attributed to the positive 
pore-water pressures escaping above the mudline due to the 
pressure gradient. 

An extreme case of pumping-erosion involved the 
Christchurch Bay Tower located in 8.4 meters of water. The 
platform was subject to severe storm conditions. Due to the 
absence of protective scour skirts the platform was 
undermined by scour, which lead to a free rocking motion of 
the structure, resulting in pumping-erosion on the 
foundation soil. 

Pumping-erosion of gravity structures is now rare, due 
to their size and the use of scour skirts. Most settlement 



of gravity structures, as described by DiBiagio (2), is 
attributed to consolidation settlement. 

Teramoto et al . (13) discussed the scour encountered with 
various sit-on-bottom type of structures. The study 
investigates the various scour patterns and contours for 
various footing designs and configurations. 

However, pumping erosion is more applicable to 
lightweight mobile structures (i.e. oceanographic research 
instrumentation packages which can be deployed from research 
vessels) and is the focus of this topic. 

Several works have dealt with wave induced pore water 
pressures. Approaches have varied from considering the soil 
skeleton to be rigid, to being compressible and pore water 
as being compressible or incompressible. Governing 
equations were usually derived from Darcy's Law of flow 
through a permeable bed or Biot's Law for three dimensional 
consolidation for a poro-elastic material. In addition, 
various studies have evaluated whether residual pore-water 
pressures exist in cohesionless soils or whether the pore 
pressures are merely transient. 

Old inziel and Brink (9), indicate that an upward flow 
of water through a porous sand bed reduces the apparent 
weight of the sand particles and therefore reduces the sand 



particles' stability. The results of the s.tudy conclude 
that upward flow (blowing) through the bed increases the 
rate of sediment transport. The study relates the upward 
blowing forces to pressure gradients in the bed that depend 
on seepage velocities as defined by Darcy's Law. 

Martin and Aral (5), 1971, indicates that the seepage 
force on a surface grain is only 50 1 of that for an embedded 
grain. However, the study concludes that it is clear that 
upward seepage will reduce the stability of surface grains, 
whereas the downward seepage will increase stability. 

Martin (4) in 1974, demonstrated that bed failure could also 
result from a horizontal pressure gradient. 



CHAPTER III 



PROCEDURES AND INSTRUMENTATION 



When non-fixed seafloor systems are subjected to wave 
and current forces of sufficient levels, the structure will 
oscillate. The wave pressure forces, in conjunction with 
the current forces will lift the structure and then will 
load it again on the seafloor. The response of the 
structure and subsequent subsidence are a function of 
several geometric, hydrodynamic and kinematic factors. 

Due to the amount of limited data concerning pumping- 
erosion subsidence a simple plate type footing employing 
purely vertical oscillations was chosen. The tests were 
conducted in a steel tank with a glass observation window. 
The dimensions of the tank were 8 feet long, 4 feet wide and 
4 feet deep. The glass observation window was 57 inches 
long and 3/4 inches thick. A sand bed, 11.5 inches deep was 
placed in the steel tank. The test section (4 feet by 4 
feet) was 11.5 inches deep. The remainder of the tank had 
5.5 inches of sand overlying a false bottom (Figure 2). 

The sand utilized for the test was an Ottawa sand with 
a specific gravity of 2.65 and a mean grain diameter (d50) 
of ,22mm (0.087 inches) (Figure 3). The sand had a static 
(pluivated) submerged angle of repose of 16. 5o. 
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Sediment Test Tank 
w/False Bottoms 

FIGURE 2 
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FIGURE 3 



Two steel plates were utilized in the experiments. A 
12 inch square plate (1 square foot), 3/4 inch, 25 pound 
plate was utilized for experiments 1 through 6. A 10 gauge 
(0.132 inch) thick, 16.97 inch square (2 square feet) plate 
attached to the 1 square foot plate was utilized for the 
experiments 7 through 12. 

I 

The plate was uplifted by a cable and pulley system 
coupled to a 1/2 horsepower D.C. variable speed motor via a 
rotating cam follower that engaged lever arm attached to the 
cable (Figure 4). This particular arrangement allowed the 
plate to be uplifted from the sandbed to a predetermined 
amplitude and then the plate was released and allowed to 
free fall onto the test bed. The next cycle would engage 
the cam follower on the lever arm lifting the plate again, 
repeating the cycle. The plate was thus able to follow the 
contours of the scour caused by the pumping-erosion. 

The amplitude of oscillation (1 or 0.8 inches) was 
obtained by eccentrically locating the cam follower from the 
center of rotation of the motor shaft. The selection of the 
amplitudes was restricted by equipment limitations. 

Oscillation periods of 5, 10 and 20 seconds were 
selected to simulate general storm wave conditions that 
might be experienced by seafloor systems. 



Two Block Pulley 







CL) 

4-J 


l , » 




03 


; ‘ • 


— 


\ 


> 

3 

i A 


9 1 




111 » , 




1 

CH C- y 

y— n i / 


/ V \ 


CD <d 

-n 


UJ ( 

03 L> 

c — « 




w JLJ 

•<- =3 

— * 1 


1— T3 


\ '' ' 


j f - ~ 

CD 


T3 xT^ 
C CD 


i ; ; ■ 

J ; 


t> 


<T3 CQ 
UO 







15 



Experimental Set-Up 

FIGURE /| 



Pore-water pressures in the test bed were measured by a 
Statham P-22S-350 bi-directional pressure transducer (2 
psid). The transducer was placed outside of the tank at the 
same equivalent hydrostatic depth as the sand bed. The 
transducer detected pressure fluctuations through the active 
port via a 1/8 inch clear PVC tubing coupled to a needle 
shaped copper probe fitted with a # 200 mesh screen to 
restrict the entrance of sand particles into the transducer 
port. The transducer probe was centered under the plate- 
footing at a depth of 0.276 inches (7mm). 

To insure proper measurements, the active transducer 
port, PVC tubing, brass coupling and sensing probe were 
thoroughly cleaned with a strong detergent to remove any 
dirt, grease or oil that may have accumulated. After 
rinsing, the active pressure port and extensions were 
carefully flushed and filled with de-aired water. This 
procedure was followed to eliminate the possibility of 
entrapped air within the pressure sensing system. 

A static load test was utilized to calibrate the 
pressure transducer. Readings were taken with the 
transducer submerged in still water at a predetermined 
depth. At a second predetermined depth readings were taken 
again. Since the pressure relationship is a linear, an 
average calibration factor of 0.043 psi/mv was obtained. 



The 5 volt D.C. excitation voltage was supplied by a 
HP17403A pre-carrier amplifier. The output was recorded by 
a HP7402A dual channel strip chart recorder. 

To record the near bed shear velocities at the plate 
edge, a Thermo-Systems Inc (TSI), hot film anemometer system 
was employed. The voltage output of the anemometer is 
proportional to heat loss across the cylindrical filament 
(film) as a result of the fluid flow across the filament. 

The actual velocities were determined from calibration 
curves developed prior to each experiment. 

The calibration was performed by attaching the hot film 
probe to a motorized trolley carriage which rode on a rail 
system above a 120 foot wave flume. The carriage was 
operated at several different speeds over a predetermined 
distance. The voltage output from the TSI anemometer 
monitor was displayed on a Beckman 200 digital voltmeter and 
recorded. A calibration curve of millivolts versus ft/sec 
was thus developed. Temperature corrections based on 
Reynolds Number similarities were employed to compensate for 
the difference in water temperatures between the wave flume 
and the test tank. 
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The upward plate velocities (VPU) and the downward 
plate velocities (VPD) were determined from the hot-film 
anemometer velocity record. 

The scour depth readings were recorded manually by 
observing a marker attached to a stanchion affixed to the 
plate traveling across a stationary rule. Horizontal 
deflections were minimized by extending a guide tube over 
the square tubing shaft welded to the plate surface. 

Twelve experiments were conducted at various 
combinations of period (T), amplitude of oscillation (A), 
and plate dimensions (area (AP); weight (W) ; and thickness 
(T) ). 



CHAPTER IV 



EXPERIMENTAL RESULTS 

A total of 14 experiments were conducted. Two 
preliminary experiments were performed to observe the 
effects of the transducer probe position under the plate. 

As one might expect, the highest negative pore-water 
pressures were recorded at the center of the plate. 
Therefore, the center of the plate was chosen for the probe 
placement. It was the goal of these experiments to obtain 
a general understanding of the mechanisms involved in the 
phenomenon of pumping-erosion. The experiments were also 
designed to investigate whether negative pore-water 
pressures and possible fluidization of the foundation bed 
played a significant role in pumping-erosion. Therefore, it 
was necessary to discern the effects of the experimental 
parameters in pumping-erosion subsidence. Following is a 
table listing the controlled variables for each test run: 



TABLE 1 - Controlled Parameters 



EXP 

NO. 


PERIOD 

(T) 

sec • 


AMPLITUDE 

(A) 

inches 


PLATE 

AREA 

(AP) 

ft2 


PLATE 

WEIGHT 

(W) 

lbs 


PLATE 

THICKNESS 

(TH) 

inches 


1 


10 


1 


1 


25 


.75 


2 


20 


1 


1 


25 


.75 


3 


5 


1 


1 


25 


.75 


4. 


20 


1 


1 


25 


.75 


5. 


10 


1 


1 


25 


.75 


6. 


5 


1 


1 


25 


.75 


7. 


..20 


2 


.80 


36 


.132 


8. 


10 


2 


.80 


36 


.132 


9. 


5 


2 


.80 


36 


.132 


10. 


20 


2 


.80 


36 


.132 


11. 


10 


2 


.80 


36 


.132 


12. 


5 


2 


.80 


36 


.132 



The soil used was an Ottawa sand with a mean grain 
diameter of .22mm. The average bed density of the 
foundation bed (the area directly under the plate) was 3.30 
slugs/ft3. In addition the controlled parameters listed in 
Table 1, the following variables were also recorded: number 

of cycles(CYCS) , positive pore water pressures (PPWP), 
negative pore water pressures (PPWS), plate exit velocities 
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(VE), plate upward velocities (VPU), and plate downward 
velocities (VPD). 

Figures 1 through 12 in Appendix A plot the subsidence 
versus time and cycles. Figures 13 through 24 plot the 
subsidence versus the log of time and cycles. 

Figures 25, 26 and 27 are group plots of subsidence 
versus cycles and Figures 28, 29 and 30 are group plots of 
subsidence versus time for the periods of 5, 10 and 20 
seconds respectively. Figure 31, plots the subsidence 
curves versus cycles (Figure 32 versus time) for the first 
six experiments (1 sf plate). Figure 33, plots the 
subsidence versus cycles (Figure 34 versus time) for the 
last six experiments (2 sf plate). 

Figures 5A through 5L show typical strip chart output 
records for the 12 experiments. The pore water pressure 
r ecord output is recorded on channel 1 and the hot film 
anemometer output record on channel 2. 

As noted in the experimental record, in Appendix B, the 
output of the hot film anemometer became unreliable as the 
vater temperature increased above 86oF. Thus VE was not 
:onsidered in the analysis described in Chapter V. 
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PRESSURE AND VELOCITY RECORD EXPERIMENT 1 




PRESSURE AND VELOCITY RECORD EXPERIMENT 2 




PRESSURE AND VELOCITY RECORD EXPERIMENT 3 




VELOCITY 

PRESSURE AND VELOCITY RECORD EXPERIMENT 4 
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PRESSURE AND VELOCITY RECORD EXPERIMENT 5 




en 



o 

o 

HP 

w 

> 



27 



PRESSURE AND VELOCITY RECORD EXPERIMENT 6 
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PRESSURE AND VELOCITY RECORD EXPERIMENT 7 
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PRESSURE AND VELOCITY RECORD EXPERIMENT 8 




PRESSURE AND VELOCITY RECORD EXPERIMENT 9 
FIGURE 5?I 
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PRESSURE AND VELOCITY RECORD EXPERIMENT 10 




PRESSURE AND VELOCITY RECORD EXPERIMENT 11 



CHAPTER V 



DISCUSSION OF RESULTS 



As one reviews the results depicted by the figures 
(Appendix A), certain trends appear. The first observation 
is that the experiments that produce the most rapid 
subsidence have the larger plate areas, greater weights, 
larger amplitudes of oscillation, and shorter periods. 
Secondly, the experiments that subside more rapidly also had 
higher negative pore water pressures and thus higher 
liquefaction indexes as derived by Sleath (12), where the 
point of fluidization is regarded as : 

_L 6ft °1 ~ p 

pg <$y p 

p = water density 
p = bed density 

<Sp = pore-water pressure through the transducer depth dy 

<$y 



Failure of the bed could also result from too large a 
horizontal gradient : 

1 6£_ = p t --o tan 0 

“T-g" <Sy p 

where e = internal angle of friction 

As the plate is uplifted, as depicted in Figure 6, 
negative pore water pressures are developed locally in the 
foundation soil directly under the plate. The data indicate 
that the bed is indeed fluidized. Assuming that the pore- 
water pressure is incompressible and that the grain skeleton 
is compressible and deformable, the action of local 
fluidization of the bed would than cause distortion of the 
upper grain layers (see Figure 7). 

Some of the loose grains were observed to be carried 
into suspension by the turbulent flow of fluid rushing in 
under the plate in response to the local negative pressures. 
The strip chart record reflects the velocity response to the 
pressure fluctuations which are in agreement with the 
Bernoulli equation : 

p = -j pV2 - gx + C 
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Upward Plate Cycle 
' FIGURE 6 
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Distorted Grain Skeleton 
FIGURE 7 



During the downward cycle of the plate (Figure 8), the 
loosened grains as well as the sediment in suspension (as a 
result of the uplifting of the plate) can then be easily 
transported (Figure 8) at a much lower shear velocities (u*) 
then would be required for normal initiation as predicted by 
the Shields' diagram (6). An interesting observation, 
derived from the pressure and velocity records, is that the 
slower periods have more pressure fluctuations and thus more 
velocity fluctuations but with lower negative pore water 
pressures than do the shorter period cycles which have 
greater negative pore water pressures but with a faster 
change in pressure and velocity with fewer fluctuations. 

By careful observation three interesting details were 
discovered. First, a small ridge of sand developed along 
and under the perimeter edge of the plate. The second 
observation is that a trough was developed immediately 
adjacent to the perimeter edge of the plate (Figure 9). 

Thirdly, slope instability and sloughing was observed during 
the upstroke of the plate. These seemingly minor details 
may play a significant role in the pumping erosion process. 

The small ridge under the perimeter of the plate 
probably results from the outflow of water and transported 
material attempting to exit from the underside of the plate 
during the last fraction of a second before plate impact. 

This small ridge of sand may account for the large negative 
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Bed Profile Details 



pore water spikes at the end of each cycle (Figure 10) in 
one of two ways. First, if the foundation bed is considered 
to be elastic, the negative pore water pressure may be due 
to plate-soil bounce with the small ridge acting to seal the 
plate so that the pressure response is localized. The 
second possibility, is that as the plate impacts on the 
foundation soil, the plate undergoes slight flexure and 
causes the negative pressure excursion. Again, the sand 
ridge would tend to seal off the plate to cause the event to 
be a local effect. 

The trough is caused by the jetting of water from under 
the plate (Figure 9). During the uplifting of the plate, 
the lower portion of the slope became unstable and sloughed 
onto the aforementioned trough. The entire slope then 
cascaded down to form a new equilibrium slope. The slope 
failure was probably due to the pressure/velocity field near 
the edge of the plate. As the water rushes in, responding 
to the negative pressures under the plate, sand grains at 
the base of the slope could be enveloped by the turbulent 
flow causing a slope failure. An alternate explanation, is 
that there is flow out of the slope in response to the 
pressure fluctuation, which would also cause the slope 
failure. 
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TYPICAL PORE- WATER PRESSURE RECORD 
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NEGATIVE PORE-WATER PRESSURE SPIKE 



In any event, the sand from the slope failure, which 
should be large grains due to preferential transportation 



(winnowing), loosely fall into the trough only to be 
expelled and transported during the downward stroke of the 
plate. 

Topographical plots of the sand contours for the 12 
experiments are in Appendix C. Of particular interest is 
the dynamic angle of repose of the scour mound slopes as 
1 isted in Table 2. 

TABLE 2 DYNAMIC ANGLES OF REPOSE 
EXP. NO. DYNAMIC ANGLE OF REPOSE 



1 


21 


2 


24.3 


3 


24.9 


4 


24.2 


5 


23.4 


6 


25.8 


7 


23.3 


8 


22.2 


9 


23.8 


10 


22.2 


11 


22.4 


12 


22.7 



Figure 11, ranks the experiments by the average initial 
liquefaction index. The ranking is in accordance with the 
observed subsidence rates (i.e. the tests with more rapid 
subsidence have higher liquefaction indexes). A 
relationship between the period, amplitude, plate weight and 
plate area ( A2Ap/WT2 ) was then developed and ranked 
(Figure 12). Curiously enough, the ranking was the same as 
the liquefaction index except for the positions of 
experiments 2 and 4, which were reversed. 

Although some general trends were developed from the 
plotted data and figures, a more sophisticated numerical 
analysis is required to ascertain the relationships between 
the experimental parameters. 

MULTIVARIATE PROCEDURES 



To determine the relationship and relative importance 
of the variables the Statistical Analysis System, or SAS, 
was employed. The SAS techniques were used to analyze the 
data because of the statistical integrity and ease of use. 
All of the statistical methods employed are explained in the 
SAS User's Manual : Statistics (11). The data from all 12 
experiments were compiled into 234 observations, using 15 
variables. The dependent variable was subsidence (SU8SID). 
The 14 independent variables were : period (T), amplitude of 
oscillation (A), plate area (AP), plate weight (W), plate 
thickness (TH), bed density (8RH0) , kinematic viscosity 
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(VISC), number of cycles (CYCS), negative pore water 
pressures (PPWS), liquefaction index (LI), positive pore- 
water pressures (PPWP), plate upward velocity (VPU), plate 
downward velocity (VPD) and plate exit velocity (VE). Table 
3 summarizes the range of values for each of the variables. 

A list of all the variable values can be found in Appendix 
B. 



The SAS procedure, PROC CORR, was first used to 
determine simple linear correlations between the variables 
in the entire data record and are listed in Table 4. Due to 
the high correlation between the area of the plate (AP) and 
the plate weight (W) and plate thickness (TH), W and TH were 
deleted form further analysis. Also LI was deleted from 
analysis due to the high cross correlation with PPWS. 
Furthermore, PPWP was not considered for further analysis 
because of its low correlation (-0.016) with subsidence. 

This is an interesting event as data suggests that positive 
pore-water pressures have little to do with subsidence. 

The next step was to determine the strength of the 
linear relationship determined by computing R2 , the 
square of the multiple correlation coefficient (R-SQUARE in 
SAS). R2 also known as the coefficient of determination, is 
interpreted as the portion of the variability that has been 
accounted for by the regression formula. The closer R2 was 
to 1, the better the model equation fit the data. However, 



SUBSID 234 0.87252137 0.49662039 204.17000000 
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TABLE 4 (cont. 



a large R2 did not necessarily mean that the proper model 
was chosen for the data. Therefore, R2 should be used with 
caution, since it is always possible to make R2 large by the 
addition of more variables to the model. The R-SQUARE 
routine was configured to give R2 and C(p) values for all 
possible combinations of 3, 4 and 5 variables. The results 
are listed in Appendix D. C(p), Mallow's statistic, is a 
criteron related to the mean square of a fitted value (7). 

As C(p) attempts to indicates the bias of the prediction, a 
small C(p) value is generally sought. The "best fit" model 
is determined by when C(p) first approaches p on a plot of 
C(p) versus p. 

Next, the remaining 10 variables were incorporated into 
a SAS routine called STEPWISE. STEPWISE is a model building 
program which chooses variables that are to be included into 
the model. Ideally, a variable selection process should 
provide for any given number of independent variables, a 
subset whose equation of predicted values produce the 
minimal residual sum of the squares. In PROC STEPWISE 
variables are added or deleted one at a time, depending on 
the Forward or Backward option selection, until all 
coefficients remaining in the model are statistically 
significant (P < 0.0001). 



The STEPWISE procedure, using the backward elimination 
method, gave results that are suspect. It indicated that 
the best model was produced when all the variables were 
significant at the P < 0.001 level when the variable AP 
(plate area) was removed. This result is contrad ictory to 
the forward selection process which adds variables in the 
order of there significance. The forward selection process 
ranked the variables in order of significance in the 
following manner: 

TABLE 5 FORWARD SELECTION RANKING 



Variables 


C(P) 


F 


PR0B>F 


R2 


CYCS 


185.049 


381.754 


0.0001 


0.6220 


VPD 


28.146 


143.304 


0.0001 


0.7667 


vise 


20.663 


8.842 


0.0033 


0.7754 


AP 


16.634 


5.737 


0.0174 


0.7808 


T 


14.272 


4.209 


0.0414 


0.7848 


PPWS 


11.090 


5.091 


0.0250 


0.7895 


A 


9.307 


3.761 


0.0537 


0.7930 


VPU 


9.000 


2.307 


0.1302 


0.7951 



The significance of the ranking is apparent in the 
manner in which the forward selection was constructed, in 
that it added the variables in the following order : CYCS, 

VPD, VISC, AP, T, PPWS, A and VPU. Each additional variable 
improved the overall fit of the model. 

The strength of a linear relationship, an indication of 
how well the model fits the data, is indicated by the R2 
value. The F-value is used by SAS to test the null 
hypothesis, which tests the linear association between Y and 
XI, X2, . . . Xk . If the F-value is not significant (PROB F 
= 0.0001), for a system on "n" degrees of freedom (DF), the 
null hypothesis is accepted indicating that the model does 
not explain a significant portion of the variability. That 
is to say, as the F-value increases a greater portion of the 
variability is explained by the model (i.e. the spread of 
the model is lessened). 

The PROB >F indicates the significance of the model 
(or the strength of the variables as in the R2 output). The 
lower the value of PR0B>F , the more significant the model . 
A p-value of 0.0001 is highly significant. 

The strength of the individual variables in indicated 
by, T and PR0B>T, which denote the same inference as F and 
PR0B>F. 



Based on the above results, PROC REG was performed on 
the 8 remaining independent variables. Variables A 
(amplitude) and V PU (plate upward velocity) were then 
removed and PROC REG was run again. The results are shown 
in Tables 6 and 7. 

Two variables, FI and F2 were then created from the 
combination of the 8 variables where : 



FI = {(144 PPWS VPU AP/( g2T))}l/3 

F2 = { ( 1 / 2 VPD 2 A W ) / (144 VISC2 g2 T2 )}l/2 

F3 and F4 were created from the log of FI and F2 taking 
into account the number of cycles. PROC CORR, PROC RSQUARE 
and PROC REG were then utilized on FI, F2 F3, and F4. The 
results are listed in Tables 8, 9 and 10. 

Plots of the models (Figures 13, 14 & 15) listed in 
Tables 6, 7 and 10, graph the predicted values from the 
models against the actual values. 
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Multivariate Results 



The SAS multiple regression was used to determine if 
the dependent variable Y (subsidence) was linearly related 
to two or more of the independent variables XI, X2,...Xk . 
The solution takes the form of : 

Y = a + blXl + b2 X 2 +...+ bkXk 

where "a" denotes the intercept term when X's equal 
zero, bl , b2,...bk are the partial regression coefficients 
of XI, X2,...Xk , respectively. 

Below is listed a summary of MODI, the 8 variable 
model; M0D2, the 6 variable model; and M0D3, the 2 variable 
model containing FI and F4. 

TABLE 9 MODEL SUMMARY 

SUM OF 

MOD # VARS R2 F-VALUE SQ. RESIDUALS 



1 


8 


0.7951 


109 


11.78 


2 


6 


0.7895 


142 


12.09 


3 


2* 


0.6618 


226 


19.43 



* 2 terms containing 8 variables 



From Table 11, the 6 variable model was selected as the 



best model where: 

Y = -1.72105 - 0.11987 T + 0.7103 AP + 0.00536 CYCS + 
1.85197 PPWS + 21645.59 VISC + 0.16325 VPD 

The determination of the "best" linear model is not as 
important to this research as is the significance of the 
variables. From Tables 5 and 7, it can be seen that the two 
most significant variables are the number of cycles (CYCS) 

(F = 381.754, T = 18.39, p-value = 0.0001) and downward 
plate velocity (VPD) (F = 143.304, T = 10.225, p-value 
=0.0001). The remaining variables listed in their order of 
significance are : viscosity (VISC), plate area (AP), 
Negative pore water pressure (PPWS) and period (T). 



CHAPTER VI 



SUMMARY AND CONCLUSIONS 



The main objective of this research was to determine 
which mechanisms were involved in the phenomenon of pumping 
erosion. The secondary objectives, involved the evaluation 
of the relative importance of each of the 14 independent 
variables, and in particular, the significance of negative 
pore-water pressures. 

Of the original 14 independent variables, 6 variables 
(cycles (CYCS) , downward plate velocity (VPD), viscosity 
(VISC), plate area (AP) , negative pore water pressure 
(PPWS), and period (T)) were considered significant in the 
model constructed utilizing SAS. The number of cycles 
(CYCS) and the downward velocity of the plate (VPD) were 
highly significant and is in agreement with the accepted 
understanding of pumping-erosion and scour. 

The result of negative pore water pressure (PPWS) 
remaining significant throughout the analysis supports the 
hypothesis that negative pore water pressures are a 
contributing factor in pumping-erosion. This is 
understandable if the experimental results of the 



Liquefaction Index (fluidization index) associated with the 
negative pore water pressures is accepted. 



This leads to the conclusion that the fluidization of 
the foundation soil, due to the negative pore water 
pressures developed during the upward cycle of the plate- 
footing is a precursor to the scour developed during the 
downward cycle of of the plate- footing and significantly 
contributes to the overall pumping-erosion subsidence . 

As with many initial investigations, this research 
poses more questions than it answers. Such as : 

1) If negative pore-water pressures and the 
subsequent fluidization are significant to the 
phenomenon of pumping-erosion, then what is the 
distribution of the negative pore water 
pressures under the footing ? 

2) What is the distribution of velocities under 
the plate footing ? 

3) Can the footing be designed to alter the 
pressure and flow nets under the footing to 
reduce subsidence ? 



4) What effect does the shape of the scour mound 
have on the negative pore water pressures ? 

5) How would the results vary in other 
cohesionless and cohesive soils ? 

6) What are the effects on pivotal oscillations 
as opposed to vertical oscillation ? 



RECOMMENDATIONS 



To obtain a better understanding of pumping-erosion, 
more research is needed in this area. If the results are 
valid, the effect of negative pore water pressures may 
affect all dynamic structure-mar ine soil interactions 
including submarine pipelines. 

For purposes of further research the following 
recommendations are proposed : 



1) Employ hydraulic or pneumatic rams to load the 
soil at a constant rate with a constant 
amp! itude. 
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2) Develop a method to record the dynamic 
formation of the scour mound slope. 



3) Conduct the experiments with different grain 
sizes of cohesionless soils. 



4) Conduct the experiments with cohesive soils. 
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APPENDIX A 



SUBSIDENCE CURVES 



UBSlDENCE (Inches) SUBSIDENCE (Inches) 



Subsidence vs Time 

EXPERIMENT 2 Ts 20 sees A= I inch 




Subsidence vs Cycles 



EXPERIMENT 2 7s 20 sees A- l inch 




FIGURE 2 



SUBSIDENCE (Inches) SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 3 T=. 5 secs A= 1 inch 




Subsidence vs Cycles 

EXPERIMENT 3 T= 5 secs A= 1 inch 




FIGURE 3 



SUBSIDENCE (Inchas) SUBSIDENCE (Inchaa) 



Subsidence vs Time 

EXPERIMENT 11 T= 10 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 

EXPERIMENT 11 T= 10 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




FIGURE 11 



SUBSIDENCE (inches) SUBSIDENCE (Inches) 



Subsidence vs Time 

EXPERIMENT 12 T= 5 s«cs A= 0.8 inch Ap= 2sf Wp = 36lbs 




Subsidence vs Cycles 



EXPERIMENT 12 T= 5 sees As 0.8 inch Ap= 2s 1 Wp= 36lbs 




SUBSIDENCE (Inches) SUBSIDENCE (inches) 



Subsidence vs Time 

EXPERIMENT 1 T= 1 0 secs A= I inch 




Subsidence vs Cycles 

EXPERIMENT 1 T= 1 0 secs A- 1 inch 




FIGURE 1 



SUBSIDENCE (inches) SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 1 T= 1 0 secs A= I inch 




Subsidence vs Cycles 

EXPERIMENT 1 T= 10 secs A= 1 inch 




FIGURE 13 



SUBSIDENCE (inches) SUBSIDENCE (Inches) 



Subsidence vs Time 

EXPERIMENT 2 T= 20 secs A= I inch 




Subsidence vs Cycles 

EXPERIMENT 2 T= 20 secs A= I inch 




FIGURE 14 



SUBSIDENCE (Inches) SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 3 T= 5 sees A= 1 inch 




Subsidence vs Cycles 



EXPERIMENT 3 T= 5 secs A= 1 inch 




FIGURE 15 



SUBSIDENCE (inches) SUBSIDENCE (Inches) 



Subsidence vs Time 



EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= 1 sf Wp= 25lbs 




Subsidence vs Cycles 

EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= 1 sf Wp= 25lbs 




FIGURE 16 



SUBSIDENCE (Inch**) SUBSIDENCE (Inches) 



Subsidence vs Time 



EXPERIMENT 5 T= 10 secs A= 0.8 inch Ap= Isf Wp= 25lbs 




Subsidence vs Cycles 

EXPERIMENT 5 T= 10 secs A= 0.8 inch Ap= Isf Wp= 25lbs 




FIGURE 17 



SUBSIDENCE (Inches) SUBSIDENCE (inches) 



Subsidence vs Time 

EXPERIMENT 6 T= 5 secs A = 0.8 inch Ap= Isf Wp= 25lbs 




Subsidence vs Cycle 

EXPERIMENT 6 T= 5 secs A= 0.8 inch Ap= Isf Wp= 25lbs 




figure is 



SUBSIDENCE (inches) SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 7 T= 20 sees A= 1.0 inch Ap= 2ti Wp= 36lba 




Subsidence vs Cycles 



EXPERIMENT 7 T=s 20 sees A= 1.0 inch Ap= 2sf Wp= 36lba 




FIGURE 19 



SUBSIDENCE (Inches) SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 8 T= 10 s«cs A* 1.0 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 

EXPERIMENT 8 T= 10 s«es A= 1.0 inch Ap= 2sf Wp= 38lbs 




FIGURE 20 



SUBSIDENCE (lnch«») SUBSIDENCE (lnch«») 



Subsidence vs Time 

EXPERIMENT 9 T= 5 secs A= 1.0 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 

EXPERIMENT 9 T= 5 secs A= 1.0 inch Ap= 2sf Wp= 36lbs 




FIGURE 21 



SUBSIDENCE (Inches) SUBSIDENCE (Inches) 



Subsidence vs Time 



EXPERIMENT 10 T= 20 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 

EXPERIMENT 10 T= 20 secs A= 0.8 inch Ap= 2s 1 Wp= 36lbs 




FIGURE 22 



SUBSIDENCE (lnch«») SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 11 T= 10 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 



EXPERIMENT 11 T= 10 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




FIGURE 23 



SUBSIDENCE (inches) SUBSIDENCE (inches) 



Subsidence vs Time 

EXPERIMENT 12 T= 5 sees A= 0.8 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 

EXPERIMENT 12 T= 5 secs A= 0.8 inch Ap= 2sf Wp= 36lbs 




FIGURE 24 
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STEP 1 VARIABLE CYCS ENTERED R SQUARE = 0.62199808 C(P) = 185.04939267 
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BOUNDS ON CONDITION NUMBER: 1.140449. 9.860305 



STEP 4 VARIABLE AP ENTERED R SQUARE = 0.78084502 C(P) * 16.63407610 
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STEP O ALL VARIABLES ENTERED R SQUARE = 0.79508359 C(P) = 9.00000000 
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MAXIMUM R-SOUARE IMPROVEMENT FOR OEPENDENT VARIABLE SUBSIO 
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SUBSIDENCE (lnch«») SUBSIDENCE (Inches) 



Subsidence vs Time 



EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= Isf Wp= 25lbs 




Subsidence vs Cycles 



EXPERIMENT 4 T= 20 secs A= 0.8 inch Ap= Isf Wp= 25lbs 




FIGURE 4 
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SUBSIDENCE (Inches) SUBSIDENCE (Inches) 



Subsidence vs Time 

EXPERIMENT 5 T= 10 secs A= 0.8 inch Ap= 1 sf Wp= 25lbs 




Subsidence vs Cycles 



EXPERIMENT 5 T* 10 secs A= 0.8 inch Ap= 1 sf W p= 25lbs 




SUBSIDENCE (lnch«») SUBSIDENCE (inches) 



Subsidence vs Time 



EXPERIMENT 6 T= 5 secs A = 0.8 inch Ap= Isf Wp= 25lba 




Subsidence vs Cycle 



EXPERIMENT 6 T= 5 secs A= 0.8 inch Ap= 1 sf Wp= 25lbs 




FIGURE 6 
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Subsidence vs Time 



EXPERIMENT 7 T= 20 sees A= 1.0 inch Ap= 2s i Wp= 36lbs 




Subsidence vs Cycles 



EXPERIMENT 7 T= 20 secs A= 1.0 inch Ap= 2sf Wp= 36lbs 




FIGURE 7 
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Subsidence vs Time 



EXPERIMENT 8 T= 10 s«cs A= 1.0 inch Ap= 2s i Wp= 36lbs 




Subsidence vs Cycles 



EXPERIMENT 8 T= 10 s«cs A= 1.0 inch Ap= 2sf Wp= 36lbs 




FIGURE 8 
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Subsidence vs Time 



EXPERIMENT 9 T= 5 secs A= 1.0 inch Ap= 2sf Wp= 36lbs 




Subsidence vs Cycles 



EXPERIMENT 9 T= 5 secs A= 1.0 inch Ap= 2sf Wp= 36lbs 




FIGURE 9 
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